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We report on a search for charged massive long-lived particles (CMLLPs), based on 5.2 fb^' of integrated 
luminosity collected with the DO detector at the Fermilab Tevatron pp collider. We search for events in which 
one or more particles are reconstructed as muons but have speed and ionization energy loss (dE/ dx) inconsistent 
with muons produced in beam collisions. CMLLPs are predicted in several theories of physics beyond the 
standard model. We exclude pair-produced long-lived gaugino-like charginos below 267 GeV and higgsino-like 
charginos below 217 GeV at 95% C.L., as well as long-lived scalar top quarks with mass below 285 GeV. 

PACS numbers: 13.85.Rm,14.80.Ly 



We report on a search for massive particles, which are 
electrically charged and have a lifetime long enough to es- 
cape the DO detector before decaying. Charged massive long- 
lived particles are not present in the standard model (SM) 
nor are their distinguishing characteristics (slow speed, high 
dE/dx) relevant for most high energy physics studies. Al- 
though the distinctive signature in itself provides sufficient 
motivation for a search, some recent extensions to the SM sug- 
gest that CMLLPs exist and are not yet excluded by cosmo- 
logical limits [1, 2J. Indeed, the standard model of big bang 
nucleosynthesis (BBN) has difficulties in explaining the ob- 
served lithium production. The existence of a CMLLP that 
decays during or after the time of BBN could resolve this dis- 
agreement yl] . 

We derive cross section limits for CMLLPs and compare 
them to theories of physics beyond the SM. In most super- 
symmetric (SUSY) models the lightest SUSY particle is as- 
sumed to be stable. Some SUSY models predict that the next- 
to-lightest supersymmetric particle (NLSP) can be a CMLLP. 
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In this Letter we explore models which include a chargino 
as a NLSP whose lifetime can be long if its mass differs 
from the mass of the lightest neutralino by less than about 
150 MeV[45 JJ. This can occur in models with anomaly- 
mediated supersymmetry breaking (AMSB) or in models that 
do not have gaugino mass unification. There are two general 
cases, where the chargino is mostly a higgsino and where the 
chargino is mostly a gaugino, which we treat separately in this 
Letter. 

There are some SUSY models that predict a scalar top 
quark (stop) NLSP and a gravitino LSP These stop quarks 
hadronize into both charged and neutral mesons and baryons 
which live long enough to be CMLLP candidates [6]. Further, 
hidden valley models predict GMSB-like scenarios where the 
stop quark acts like the LSP and does not decay but hadronizes 
into charged and neutral hadrons that escape the detector 
10, Si- In general, any SUSY scenario where the stop quark 
is the lightest colored particle (which will happen in mod- 
els without mass unification and heavy gluinos) can have a 
stop CMLLP. Any colored CMLLP will undergo hadroniza- 
tion and charge exchange during nuclear interactions, which 
we discuss below. 

This search utilizes data collected between 2006 and 2010 
with the DO detector [9] at Fermilab's 1.96 TeV pp Tevatron 
Collider, and is based on 5.2 fb^' of integrated luminosity. 
We reported earlier ifioll on a similar 1.1 fb^' study, searching 
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for events with a pair of CMLLPs, each with low speed. In ad- 
dition to using the larger data sample, the present search looks 
for one or more CMLLP, rather than only for a pair, and char- 
acterizes CMLLPs with high dE / dx in addition to slow speed. 
Other searches for long-lived particles include those from the 



CDF Collaboration lU 



the CERN e+e^ Collider LEP 

ifisll , and the CERN pp Collider LHC Hfisll. 

The DO detector 1 9| includes an inner tracker with two com- 
ponents: an innermost silicon microstrip tracker (SMT) and 
a scintillating fiber detector. We find a particle's dE/dx 
from the energy losses associated with its track in the SMT. 
The tracker is embedded within a L9 T superconducting 
solenoidal magnet. Outside the solenoid is a uranium/liquid- 
argon calorimeter surrounded by a muon spectrometer, con- 
sisting of drift tube planes on either side of a L8 T iron toroid. 
There are three layers of the muon detector: the A-layer, lo- 
cated between the calorimeter and the toroid, and the B- and 
C- layers, outside the toroid. Each layer includes scintilla- 
tion counters which serve to veto cosmic rays. Thus the muon 
system provides multiple time measurements from which a 
particle's speed may be calculated. 

Because we distinguish CMLLPs solely by their speed j3 
and dE /dx, we must measure these values for each muon can- 
didate as accurately as possible. Muons from Z — > jXjX events 
studied throughout the data sample allow calibration of the 
time measurement to better than 1 ns, with resolutions be- 
tween 2-4 ns, and to maintain the mean dE/dx constant to 
within 2% over the data-taking period. From a specific muon 
scintillation counter we calculate a particle's speed from the 
time recorded and the counter's distance from the production 
point, and we compute an overall speed from the weighted av- 
erage of these individual speeds, using measured resolutions. 
The ionization loss data from the typically 8-10 individual en- 
ergy deposits in the SMT are combined using an algorithm 
that omits the largest deposit to reduce the effect of the Lan- 
dau tail and corrects for track crossing angle. We normalize 
the dE /dx measurements by requiring the dE / dx distribution 
of muons from Z ^ pLpL events to have a maximum at L Fig- 
ure [T] shows the distributions in j3 and dE / dx for data and 
background, which pass the selection criteria described below. 

The selection of a candidate CMLLP occurs in several 
steps. Because of the high pp collision rate, we employ a 
three-level trigger system to reduce the event rate to the 200 
Hz that can be recorded. The trigger system bases its decisions 
on characteristics of the event, which for the CMLLP candi- 
dates is the presence of a muon with a high momentum trans- 
verse to the beam direction {pj)- In order to reduce triggers 
on cosmic rays, there is a time window at the initial trigger 
level. This trigger gate reduces the trigger efficiency by 10% 
for CMLLPs with a mass of 300 GeV (as they will be slow 
and some will be out-of-time) and contributes significantly to 
our overall acceptance. We avoid a tighter timing gate usually 
imposed at the second level of the muon trigger by accept- 
ing an alternative requirement that the muon have a matching 
track in the SMT. 

In the standard DO event reconstruction CMLLPs would ap- 
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Figure 1 : Distributions of (a) speed jS and (b) dE /dx for data, back- 
ground, and signal (gaugino-like charginos with a mass of 100 and 
300 GeV) that pass the selection criteria. The histograms have been 
normalized to have the same numbers of events. We have adjusted 
the scale of the dE/dx measurements so that the dE/dx of muons 
from Z — 7- /i/i events peak at \. All entries exceeding the range of 
the histogram are added to the last bin. 



pear as muons, which has been verified in detail using MC 
simulations. Thus we select events with at least one well 
identified high pj muon. For a useful j3 measurement, the 
event must have scintillator hits in the A-layer and either the 
B- or C-layer; and for valid dE/dx data, we require at least 
three hits in the SMT. For an optimal tracking and momen- 
tum measurement we require the muon to be central, i.e., with 
a pseudorapidity llT6i | T] | < 1 .6 . To reject muons from me- 
son decays we impose the isolation requirement that the sum 
of the pT be less than 2.5 GeV for all other tracks in a cone 
of radius R = (A^)^ + (At])^ < 0.5. We also require that 
the total transverse calorimeter energy in an annulus of radius 
0.1 < R < 0.4 about the muon direction be less than 2.5 GeV. 
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A requirement that the z-coordinate of the muon track at the 
location of closest approach to the beam axis be < 40 cm en- 
sures that the particle passes through the SMT. 

We impose further criteria to eliminate cosmic rays. To 
select muons traveling outwards from the apparent interac- 
tion point, we require that its C-layer time be significantly 
greater than its A-layer time. We require also that the muon's 
distance-of-closest-approachto the beam line be less than 0.02 
cm. These criteria are also applied to a second muon, if any, 
in the event. In addition, for events with two muons we re- 
quire that the absolute value of the difference between each 
muon's A-layer times be less than 10 ns. To reject cosmics 
that appear as two back-to-back muons, we require for their 
pseudo-acolinearity Aa = |A0 — 27i;| > 0.05. 

Events with a muon from a W boson decay, with mismea- 
surements providing inaccurate values of the muon's j3 and 
dE I dx, constitute a potentially large background. To study se- 
lection criteria for CMLLPs, we calculate the transverse mass 
Mt Inl] and select data with Mr < 200 GeV to model the 
data in the absence of signal 11811 . We choose selection cri- 
teria that minimize the number of events surviving from this 
background sample compared to events from simulations of 
the CMLLP signal. We require that events contain at least one 
muon with pr > 60 GeV. From a separate sample of muons 
from Z ^ fifl decays we observe that the association of a 
spurious scintillator hit can result in an anomalously low )3 
value. We use an algorithm that discards such hits through 
minimizing the /d.o. f. for the j3 calculated from the dif- 
ferent scintillator layers. By comparing the effect on the back- 
ground sample with the effect on simulated signal, we choose 
to eliminate events unless the minimized speed j^/c/.o./. < 
2. Finally, we compare the muon's track direction measured 
by the muon system with that measured in the central tracker, 
and eliminate events with clearly mismatched tracks. 

To simulate signal events, we generate CMLLP candi- 
dates using PYTHIA 11191]. The long-lived stop quarks are 
hadronized using ||20|] . Because the signature of the CMLLP 
cascade decays is model dependent and difficult to simu- 
late accurately, we generate direct pair-production of the 
CMLLPs, without including cascade decays. We use the full 
DO detector GEANT L21.1 simulation to determine the detec- 
tor response for these samples (which include overlaid data- 
based pp interactions). The results are applicable to models 
with pair-produced CMLLPs with similar kinematics. 

The stop quarks are distinct since they appear in charged 
or uncharged stop hadrons, which may flip their charge as 
they pass through the detector. In the simulation approx- 
imately 60% of stop hadrons are charged following initial 
hadronization [22J, i.e., 84% of the events will have at least 
one charged stop hadron. Further, stop hadrons may flip their 
charge through nuclear interactions as they pass through ma- 
terial. We assume that stop hadrons have a probability of 2/3 
of being charged after multiple nuclear interactions and anti- 
stop hadrons, a probability of 1/2 of being charged, consistent 
with the numbers of possible stop and anti-stop hadronic final 
states 123142511 . For this analysis we require the stop hadron to 



be charged before and after passing through the calorimeter, 
i.e., to be detected both in the tracker and in the A-layer, and 
to be charged after the toroid, i.e., to be detected in the B- or 
C-layers. The probability for at least one of the stop hadrons 
being detected is then 38%, or 84% if charge flipping does not 
occur. We include these numbers as normalization factors in 
the confidence level analysis discussed below. 

Our final selection criteria is that the candidate's speed 
j3 < L Thus, we describe the background by the j3 < 1 data 
events with Mj < 200 GeV, and search for CMLLP candi- 
dates in j3 < 1 data with Mj > 200 GeV. We normalize the 
background and data samples in the )3 > 1 region, where the 
contribution of signal is negligible. Because the uncertainties 
in the speed measurements depend on the particle's Tj, due to 
detector geometry, and the distributions in rj of the muons in 
the Mt < 200 GeV sample differs from those in the Mj > 200 
GeV sample, we use the signal-free region to derive correction 
factors for the background sample that match its Tj distribution 
to that of the data. 

We utilize a boosted decision tree (BDT) li261 to discrim- 
inate signal from background. The most discriminating vari- 
ables are the CMLLP candidate's j3 and dE /dx, but we also 
include several related variables: the speed significance, de- 
fined as {I — P) /ap, the corresponding number of scintil- 
lator hits, the energy loss significance defined as {dE/dx~ 
1 )/ /(/^, and the number of SMT hits. For each mass point 
in all three signal models we train the BDT with the signal 
MC and the background, and then apply it to the data sam- 
ples. Figure [T] shows the distributions in j3 and in dE /dx for 
the data and background samples, as well as for two represen- 
tative signals. 

Systematic uncertainties are studied by applying variations 
to the background and signal samples and determining the de- 
viations in the BDT output distributions. Two of the system- 
atic uncertainties affect the shape of the BDT distribution and 
their effect is taken into account explicitly in the limit calcu- 
lation: the uncertainty due to the width of the Level 1 trigger 
gate and the uncertainty in the corrections to the MC time res- 
olution; both are applied only to signal. By examining the 
signal-like region of the BDT distributions, we find that the 
maximum (average) uncertainty is 10% (4%) for the trigger 
gate width, and 38% (7%) for the time resolution correction. 
All other systematic uncertainties affect only the normaliza- 
tion of the BDT output. The systematic uncertainties on the 
background are due to the dE /dx modeling (< 0.1%) and the 
background normalization, from the specific values used for 
the /3 (7.2%) and Mj requirements (2.2%). The systematic 
uncertainties on the signal include muon identification (2%) 
and the integrated luminosity (6.1%) fxf] . The systematic un- 
certainties associated with the corrections to the muon pj res- 
olution and to the dE / dx resolution, as well as the choice of 
PDF and factorization scale, are all below 1%. 

We obtain the 95% C.L. cross section limits from the BDT 
output distributions, constraining systematic uncertainties to 
data in background dominated regions |28ll . These limits are 
shown in Figure |2l together with the NLO theoretical signal 
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Figure 2: 95% C.L. cross-section limits as a function of mass for 
(a) gaugino-like charginos, (b) higgsino-like charginos, and (c) stop 
quarlss. The stop quark limits are displayed for the assumed charge 
flipping (charge survival probability = 38%) and for no charge flip- 
ping (charge survival probability = 84%). 



cross sections, computed with PROSPINO 112911 . The limits 
are of order 0.01 pb for directly pair-produced CMLLPs with 
masses between 200 and 300 GeV. Using the theoretical cross 
sections, we are able to exclude gaugino-like charginos be- 
low 267 GeV and higgsino-like charginos below 217 GeV. For 
stop quarks, we assume a charge survival probability of 38%, 
as discussed above, and exclude masses below 285 GeV. If 
charge flipping does not occur, we would obtain a higher mass 
limit, as indicated in Fig. |2j;. 

As seen in Figure|2l the observed limit exceeds the expected 
limit at various mass points by as much as 2.5 standard devi- 
ations, for all signals tested, due to the presence of the same 
few data events with high BDT discriminant values. This dis- 
crepancy reflects the excesses of data compared to background 
observed in Fig. [T]for the distributions both in beta (around 
0.6) and dE/dx (around 2.8). We have compared these events 
with background events and conclude that the observed ex- 
cess is consistent with a statistical fluctuation in the number 
of such events. 

In summary, we have performed a search for charged, mas- 
sive long-lived particles using 5.2 fb"' of integrated luminos- 
ity with the DO detector We find no evidence of signal and 
set 95% C.L. cross-section upper limits of order 0.01 pb for 
pair-produced CMLLPs of mass 200-300 GeV. At 95% C.L. 
we exclude pair-produced long-lived stop quarks with mass 
below 285 GeV, gaugino-like charginos below 267 GeV, and 
higgsino-like charginos below 217 GeV. These are presently 
the most restrictive limits for chargino CMLLPs, with about a 
factor of five improvement over the previous DO cross section 
Umits ifioll. 
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